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Abstract 

A uniform distribution of point defects in an otherwise perfect crystallographic structure usu- 
ally describes a unique pseudo phase of that state of a non-stoichiometric material. With off- 
stoichiometric uranium dioxide as a prototype, we show that analogous to a conventional phase 
transition, these pseudo phases also will transform from one state into another via changing the 
predominant defect species when external conditions of pressure, temperature, or chemical com- 
position are varied. This exotic transition is numerically observed along shock Hugoniots and 
isothermal compression curves in UO2 with first-principles calculations. At low temperatures, it 
leads to anomalies (or quasi-discontinuities) in thermodynamic properties and electronic structures. 
In particular, the anomaly is pronounced in both shock temperature and the specific heat at con- 
stant pressure. With increasing of the temperature, however, it transforms gradually to a smooth 
cross-over, and becomes less discernible. The underlying physical mechanism and characteristics of 
this type of transition are encoded in the Gibbs free energy, and are elucidated clearly by analyzing 
the correlation with the variation of defect populations as a function of pressure and temperature. 
The opportunities and challenges for a possible experimental observation of this phase change are 
also discussed. 
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Properties of states (or phases) of a matter and the transformation between them are one 
of the central topics of modern condensed matter physics. What we mean by saying a phase 
or a state of a matter refers to a region of space in a thermodynamic system that is chemically 
uniform, physically distinct, and (often) mechanically separable. From a microscopic point 
of view, it usually constitutes a sequence of static or dynamic distributions of particles, which 
possess specific symmetry or ordering, and the physical properties are uniform, distinct, and 
time-independent in a statistical sense. It is widely used as a classification of matters, for 
example, by differentiating them according to the state as solid, liquid, gas, or plasma, and 
subclassifying solid according to the crystallographic symmetry, labeling magnetism by the 
ordering of magnetic moments, etc. When external conditions such as pressure, temperature, 
or chemical composition are varied, a matter might change from one physical phase into 
another. This transformation is often accompanied by a sudden modification in crystalline 
structure or ordering state, and can be recognized by an abrupt change in thermodynamic 
properties. On the other hand, an abrupt change (or singularity) in derivatives of free 
energy- e.g., the specific volume, entropy, or specific heat- usually signals a reorganization 
of the particles, and thus can be used as a helpful indicator of phase transitions. Nevertheless, 
this is not always the case, especially in a complex system. In this paper, we will show that 
anomalous or quasi- discontinuous jumps take place in compression curves and other related 
thermodynamic quantities in non-stoichiometric uranium dioxide. These abrupt changes, 
however, do not correspond to any known physical phase transitions from a conventional 
point of view. Instead, we identify them as driven by a pseudo phase transition that is a 
result of sharp switch of predominant defect species. 

In principle, introduction of defects might modify or tune the physical properties of a 
material. This modification, however, is primarily a result of accumulation, and reflects the 
magnitude of defect concentrations. It is quite different from the abrupt variation in physical 
properties that arises from a symmetry or ordering change in a conventional phase transi- 
tion. Therefore, according to the conventional classification, a defective state of a material 
cannot be labeled as a distinct physical phase against the perfect one. Notwithstanding this, 
in some special situation the defect concentrations are strictly constrained by, say, chemical 
composition, and thus it has a possibility to exhibit a sudden change in predominant defect 
species when thermodynamic conditions are varied, while keeping the underlying matrix 
structure almost unchanged. In this sense, a uniform distribution of each type of defects 
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FIG. 1. (Color online) Compression curves of UO2+1E along the Hugoniot shocked from 500 K and 
GPa as a function of the effective volume per atom at different stoichiometry deviations of x=- 
0.02, -0.05, and -0.1, respectively. An isotherm of £=-0.1 at 1700 K is also plotted for comparison. 



would have some distinct physical properties, and analogous to a conventional phase tran- 
sition, might also lead to a jump in these properties when the predominant defect species is 
changed. We therefore can identify each such kind of homogeneous distribution of defects 
as a pseudo phase according to the dominant species of defects to differentiate them. As 
will be shown below, at low temperatures the variation in defect population during a change 
of pseudo phases is very sharp, and strong anomalies or quasi- discontinuities in physical 
properties due to pseudo phase transitions can be expected. 

Uranium dioxide (UO2) is an important nuclear fuel. It is also a typical nonstoichiometric 
compound, [l| in which the structure of the underlying matrix (in fluorite structure here) is 
stable over a wide range of pressure, and has a great degree of stoichiometry deviation. It 
thus provides as a good prototype to study the exotic effects of defects that manifest in 
thermodynamic properties. For this purpose, we first computed the energetics of a series of 
defective UO2 as a function of pressure to extract the contribution of defects to formation 
free energies, from which defect populations as a function of temperature and pressure were 
then calculated. After acquired these information, the Gibbs free energy as a function of 
chemical composition (or equivalently, the deviation from stoichiometry) was constructed, 
as well as the equation of state (EOS) of the non-stoichiometric system as a derived quantity. 

In particular, the total energy calculations of the defective U0 2 system were carried out 
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FIG. 2. (Color online) Anomalies in the thermal expansivity a, isothermal bulk sound velocity, 
and specific heat at constant pressure Cp along the Hugoniot that shocked from 500 K and OGPa 
in non-stoichiometric UO"2- The solid circle (at the bottom left corner) marks the experimental 
bulk sound velocity measured at ambient conditions in perfect UO2. 



with the density functional theory usin g the PAW pseudopotentialQ, 3 plus plane-wave 
scheme as implemented in VASP code. % |5[] The LSDA+ U approach was employed to treat 
the strongly correlated 5/ electronic orbitals that localized on uranium atoms, which splits 
the / band by an effective on-site Coulomb repulsion to form a Mott insulator state. [6| De- 
tails of the ab initio calculations, as well as the establishment of non-stoichiometric EOS, 
are referred to Ref. [7|. Especially, in order to eliminate electronic meta-stable states that 
introduced via strong correlations of 5/ orbitals, we adopted the quasi-annealing method 8] 
to lower the electronic energy slowly by controlling the random disturbance from the ionic 
system which acts as a kind of heat bath. Defect populations were evaluated using point de- 



fect model 



ll| and presumed a physical condition in the closed regime, which corresponds 



to the interior of a bulk material and thus no exchange of particles with the exterior can 
occur. 

With these calculations, an interesting phenomenon was observed in the hypo-stoichiometric 
regime (x < 0) of UC^+x- Figure [1] shows the Hugoniots that shocked from an initial con- 
dition of 500 K and OGPa for x = —0.02, —0.05, and —0.1, respectively, by comparing 
with an isotherm of x = —0.1 at 1700 K. A volume collapse was observed at about 39 GPa. 
Note that in calculations we fixed the underlying matrix in a fluorite structure all the time, 
therefore this volume collapse doesnot correspond to any structural phase transformation. 
Nevertheless, this exotic behavior in the compression curve is very similar to what happens 



in a first-order phase transition. 

It is worthwhile to point out that the magnitude of this volume change is proportional 
to the size of stoichiometry deviation x, which is more pronounced for x = —0.1 than that 
for x = —0.02. It is a direct consequence of accumulation effect of defect populations, 
and is quite different from a first-order phase transition in which the dependence of the 
volume collapse on stoichiometry deviation should be weak. We also observed that the sharp 
volume change becomes smooth, and even less distinguishable, at higher temperatures, as 
the isotherm at 1700 K in FigJT] illustrates, which reminds us that it is not a conventional 
phase transition. 

Furthermore, for non-stoichiometric UO2, it is not only compression curve that exhibits 
quasi-discontinuity, other thermodynamic quantities also manifest similar anomalies. Figure 
[2] plots the calculated thermal expansivity a, specific heat at constant pressure Cp, and 
isothermal bulk sound velocity that defined as pC 2 = —V (dP /dV) T along the same shock 
Hugoniots as in FigJT] Here p, P, V, and T denote the mass density, pressure, specific 
volume, and temperature, respectively. An experimental data of bulk sound velocity of 
perfect UO2 measured at ambient conditions [12] is also drawn for comparison. It can be 
seen that all of these thermodynamic quantities show a kink at the same pressure of about 
39 GPa. It is interesting to note that the anomaly in Cp is much more striking than that of 
others, and is still discernible even for x = —0.02. By comparison, the expansivity a already 
shows a nearly smooth crossover for x = —0.1, whereas at the same x the quasi-discontinuous 
change in sound velocity is just perceptible. The occurrence of quasi-discontinuity in both 
the first and second derivatives of the free energy at the transition point implies that although 
this transformation belongs to neither a first nor a second order phase transition, it has some 
features that are analogous to both of them. 

In fact, the physical mechanism of this transformation and the associated anomalies are 
completely encoded in the construction of the Gibbs free energy of a non-stoichiometric 
system, which reads [3] 



where G° is the Gibbs free energy of the perfect system and AG is that contributed by 
defects. The defect populations n are constrained by the composition equation of 




(1) 



f{n 1 ,n 2 , ...,m)=x. 



(2) 



5 



10 20 30 40 50 60 70 




Pressure (GPa) 

FIG. 3. (Color online) Variations of the shock temperature, defect populations, and the difference 
of Gibbs free energy with respect to the perfect UO2 when across the transition zone, respectively. 

To the first-order approximation, G can be written as 

G^G° + ^8G ini . (3) 

i 

Since 5Gi differs from each other, the Gibbs free energy thus could have a cross-over when 
the predominant defect switches from one type into another under a variation of external 
conditions of temperature, pressure, or chemical composition. Furthermore, since distribu- 
tion of point defects satisfies Boltzmann statistics, the population variation is in proportion 
to exp(— 8G/ksT). Therefore the crossover is smooth at high temperatures, but becomes 
sharp at low enough temperature, which leads to a discontinuity or strong anomaly in ther- 
modynamic properties. It is necessary to note that since the Gibbs free energy in Eq.Q 
is analytic, all of its derivatives must be continuous at finite temperatures in nature, even 
though some of them might appear to be quasi-discontinuous. 

Figure |3] shows the variation of shock temperature of x = —0.1, in which the kink is 
much more appreciable than in P-V curve. Also, the difference of the Gibbs free energy 
with respect to the perfect U0 2 (i.e. AG) and the defect populations as a function of shock 
pressure are plotted. It can be seen that the kink in temperature corresponds exactly to the 
point of switch of predominant defect species. At low pressure it is oxygen vacancy (0„) that 
dominates. But at a pressure higher than 39 GPa, uranium interstitial (Uj) starts to take 
over, and 0„ is suppressed completely. This switch also reflects in the Gibbs free energy as a 
function of shock pressure, and AG (i.e., that contributed by defects) displays a clear cross- 
over from one end to the another when crossing the transition zone. These demonstrate that 
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FIG. 4. (Color online) Pseudo-phase diagram of hypo-stoichiometric UO2 on T-P plane, together 
with the induced anomaly in shock temperature. Solid lines denote that for x=-0.1 and dashed lines 
for x=-0.02. Inset provides an overall view of the boundaries that separate the dominant defects 
of oxygen vacancy and uranium interstitial for x=-0.02 (open square) and -0.1 (filled square) 
respectively. 

the anomalies are results of a pseudo-phase transition in non-stoichiometric UO2 driven by 
compression, as above discussions implied. 

If we define the locus of the intersection points of defect populations (usually at the 
transition zone center) as a pseudo-phase boundary, a pseudo-phase diagram can be con- 
structed, which divides the dominant regime of different defect species. As shown in FigJU 
for hypo-stoichiometric UO2, dominates the low pressure region and U, is the main defec- 
tive component at high pressures. Also, we found that decreasing x increases the transition 
pressure slightly. The change is about several giga-pascal at high temperatures, but less 
than 1 GPa when below 500 K, as the inset of FigJU shows. Overall, the temperature de- 
pendence of the pseudo-phase boundary is weak, reflecting the fact that this transition is 
mainly driven by compression. 

It is worthwhile to point out that the width of the transition zone increases with temper- 
ature. This is a direct reflection of the dependence of defect populations on temperature. 
At a high enough temperature, the transition zone will become so wide that the cross-over 
behavior is almost indiscernible, and then the quasi- discontinuities (or kinks) in thermo- 
dynamic quantities disappear. On the other hand, it should be noted that theoretically a 
pressure- induced pseudo-phase transition at absolute zero Kelvin is possible, which, however, 
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FIG. 5. (Color online) Isothermal compression curves for x=-0.1 at temperatures of 300, 500, 700, 
1000, 1500, 2000, 2500, and 3000 K, respectively. Inset: defect population of uranium interstitial 
along the isotherms. Arrows indicate the direction of temperature increase. 



requires a large zero point motion of particles to facilitate ionic diffusion so that to make 
the switch of defect species possible. In this extreme case, the width of the transition zone 
becomes zero, and the Gibbs function has a sharp kink so that all of its derivatives become 
truly discontinuous at the transition point. It looks more like a first-order transition. The 
difference is that all derivatives of the free energy higher than the first order are divergent at 
the transition point for a first-order phase transition, while they are finite for a zero Kelvin 
pseudo-phase transition. This difference comes from the fact that for a defect type switch 
transition, the derivatives are carried out with defect populations being fixed, which is well 
defined even at the transition point. But it is not that case for a firs-order phase transi- 
tion, and only left (or right) limit can be defined. Figure |5] illustrates the characteristic of 
evolution of the crossover from being quasi-discontinuous at low temperatures to a smooth 
bridging at high temperatures. The variation of the population of U; as shown in the inset 
of Figj5]as a function of pressure and temperature, which exhibits a sharp jump at low T but 
becomes a smooth crossover at high T, explains perfectly all anomalous behavior observed 
in the Gibbs free energy and its derivatives. 

In above discussion, we did not consider the possibility of a phase segregation, or equiv- 
alently, the stability of hypo-stoichiometric UO2. Experiments indicated that the stable 



range of x at low temperature is very narrow, especially in the x < side. jl3 



Q Never- 



theless, a meta-stable U0 2+a: with x = —0.02 had been successfully synthesized at room 



temperature in laboratory. [15[ Although above discussions showed that the anomalies for 
this chemical composition is weak in the specific volume or shock temperature, a kink in 
the Cp curve is still evident, which implies that the pseudo-phase transition can be detected 
by a highly accurate measurement of Cp in UOi.gg. Furthermore, it is also possible that a 
meta-stable phase with x = —0.05 can be synthesized at high temperatures. This stoichiom- 
etry deviation would provide more noticeable anomalies, and thus is much easier to detect 
it. In addition, according to the phase diagram of UC^+x, when temperature is higher than 
1000 K, a composition of x = —0.1 might become possible. At this temperature range the 
crossover behavior should have been smeared drastically, but our calculations indicate that 
some of the thermodynamic properties, especially Cp, still exhibit a strong anomaly that 
can be detected by experiment. 

Another challenging issue for experiments to detect a pseudo phase transition is the possi- 
ble interference coming from structural transition of the underlying matrix. Our calculations 
suggested that the pseudo transition in \J0 2 + X might take place at about 39 GPa at room 
temperature. For perfect UO2, however, compression also leads to a phase transition from 
fluorite structure to cotunnite phase at ~42 GPa. [if], [l^] Considering that introduction of 
defects should weaken the stability of the matrix, the actual transition pressure of the under- 
lying matrix for off-stoichiometric compound thus might be lowered to overlap with that of 
the pseudo phase transition, and make the relevant experiments very difficult. Nonetheless, 
what we considered here is only neutral defects. Since UOo is an insulator within this pres- 

n 

sure range, the defects can be charged. This would modify the energetics of all defects. |18| 
It therefore still has the possibility that the OvUj switch occurs at a much lower pressure. 
To make a reliable assessment of this, however, requires a comprehensive evaluation of the 
energetics of charged defects in UO2, which is beyond the scope of this paper. In spite 
of that, our analyses definitely demonstrated that a pseudo-phase transition is practically 
possible in non-stoichiometric compounds. 

According to the defect-type switching mechanism, a pseudo phase transition driven by 
temperature or chemical composition rather than by compression is also possible. Similar 
anomalies can be expected for all of these transitions. In fact, the anomalies in thermal 
expansion coefficient and bulk modulus due to defect-type switch might have been observed 
experimentally in Fei^O. A large discontinuous jump in thermal expansivity and bulk mod- 



ulus were observed when the chemical composition across the point of Feo.gsO. 



19 



20| The 



underlying physical mechanism was not clearly understood at that time, and it was proposed 
might be related to defect clustering as x increased. It becomes clear now that this picture 
is perfectly in line with the pseudo phase transition just discussed above if we treat the 
defect clusters as independent species, i.e., it is the switch of defect species that leads to the 
observed quasi-discontinuous jump in thermodynamic properties. Furthermore, we noticed 
that the transition in Fei-^O is very similar to what would happen in a hyper-stoichiometric 
UO2, where the stoichiometry deviation drives a transition from point oxygen interstitial 



to COT-o cluster, as the pseudo phase diagram of Fig.3 in Ref. [2lf] implied. However, it 
should be noted that it is difficult to control the stoichiometry deviation continuously, thus 
it is still very important to study by experiment the pseudo phase transitions that driven 
by pressure, which can provide a more explicit and clear-cut physical picture. 

In summary, we presented in this paper that a pseudo-phase transition can occur in 
non-stoichiometric compounds (and alloys as well) which arises from a rapid switch of pre- 
dominant defect species driven by external conditions of pressure, chemical composition, or 
temperature, and so on. At high temperatures this transition is a smooth crossover, whereas 
it becomes quasi- discontinuous at low temperatures and has features similar to a first-order 
phase transition. A real discontinuity might be possible at K, but which requires a large 
zero point motion of ions to facilitate particle diffusion, therefore is only possible in com- 
pounds of light elements. Being analogous to a physical phase transition, anomalies also 
present in thermodynamic properties at the transition point of a pseudo phase transforma- 
tion. With hypo-stoichiometric UO2 as a prototype, we demonstrated that such kind of 
anomalies are evident in both shock Hugoniots and isothermal compression curves using ab 
initio calculations. The transition mechanism is clearly elucidated by analysis of the Gibbs 
free energy and evolution of defect populations as a function of pressure, stoichiometry de- 
viation, and temperature. We also discussed the possibility to observe this transition in 
non-stoichiometric UO2 by experiment, as well as the accompanied challenges. It is worth- 
while to point out that since electronic structure and magnetic property usually are governed 
(or tuned) by defect species, a pseudo phase transition thus indicates a rapid change of these 
properties (in addition to the thermodynamic properties that are direct derivatives of the 
Gibbs free energy). In this regard, the effect of pseudo phase transition of defects might have 
potential impacts on electronic applications, for example, in pressure-sensitive switchers that 
change the electronic state of a device from one state into another by slightly adjusting the 
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